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Abstract 

We propose a new method of creating solitons in elongated Bose-Einstein Con- 
densates (BECs) by sweeping three laser beams through the BEC. If one of the 
beams is in the first order (TEM10) Hermite-Gaussian mode, its amplitude has a 
transversal tt phase slip which can be transferred to the atoms creating a soliton. 
Using this method it is possible to circumvent the restriction set by the diffraction 
limit inherent to conventional methods such as phase imprinting. The method allows 
one to create multicomponent (vector) solitons of the dark-bright form as well as the 
dark-dark combination. In addition it is possible to create in a controllable way two 
or more dark solitons with very small velocity and close to each other for studying 
their collisional properties. 

Keywords: cold atoms, atomic Bose-Einstein condensates, solitons 
PACS: 03.75.Hh 



1 Introduction 

Atomic Bose-Einstein condensates (BECs) have received a great deal of interest since they 
were first produced a decade ago [1-3]. They can exhibit various topological excitations, 
such as vortices and solitons. The dynamics of solitons in elongated BECs [5] is the 
atom-optics version of the nonlinear propagation of light pulses in optical fibres [4]. The 
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BEC offers a remarkable freedom in terms of controlling the physical parameters such as 
dimensionality and even the sign of the strength of the atom-atom interaction [5]. 

Solitons in BECs can be of both dark and bright type. Dark solitons are formed in 
BECs with repulsive interaction between the atoms [5]. For completely dark solitons the 
condensate wavefunction is zero at the centre and changes its sign then crossing the central 
point, i.e. the condensate wave-function has an infinitely steep it phase slip at the centre [5]. 
On the other hand, the bright solitons are formed in BECs with repulsive interaction be- 
tween the atoms. The wave-function of the BEC is then localised at the centre [5] and 
goes to zero further away from this point. The dark solitons which manifests themselves as 
a density minimum moving with a constant speed against a uniform background density, 
as well as bright solitons which are shape preserving wave packets, have both been exper- 
imentally realised [6-9]. The dynamics of solitons in BECs has been extensively studied. 
This has included investigations of the stability properties [10], as well as soliton dynamics 
in inhomogeneous clouds [11], in multicomponent BECs [12,13] and in supersonic flow [14]. 
Solitons can be created in various ways with a variable degree of controllability, e.g., by 
colliding clouds of BEC [15-17] or engineering the density [18,19]. 

Traditionally dark solitons in BECs are created using phase imprinting [6,7,20-22], 
where a part of the condensate cloud is illuminated by a far detuned laser pulse in or- 
der to induce a sharp ir phase slip in the wave function. The subsequent dynamics can 
indeed develop solitons [6,7]. There are, however, some rather severe drawbacks with 
such a method of phase engineering. The resolution of the required phase slip is naturally 
restricted by the diffraction limit, i.e. the width of the phase slip should be larger than 
an optical wave-length. Furthermore the phase imprinting does not produce a density 
minimum characteristic to the dark solitons in the region of the phase change. Hence 
completely dark stationary solitons are difficult to achieve, which consequently results in 
so called grey moving solitons with a shallow density dip. 

It is of a significant interest to be able to create slowly moving, or even completely 
stationary solitons in order to test for instance their scattering properties. The shapes of 
the colliding solitons are to be preserved. In addition, a relative spatial shift is expected. 
This spatial shift, however, can only be detected for extremely slow solitons due to the 
inherent logarithmic dependence of the spatial shift on the relative velocity between the 
solitons [23,24]. The standard phase imprinting also inevitably creates phonons in the 
trapped cloud because the constructed initial state is not the exact soliton solution largely 
due to the missing density notch [6,7]. 

In this paper we show how states which have the required phase slip and density profile 
for solitons can be created by sweeping three laser beams through an elongated BEC as 
shown in Fig. [TJ If one of the beams is in the first order (TEM10) Hermite-Gaussian 
mode, its amplitude has a transversal 7r phase slip which will be transferred to atoms thus 
producing a soliton. More importantly, with a sequence of three laser beams it is possible 
to circumvent the restriction set by the diffraction limit. The laser fields reshape an atomic 
wave-function so that it acquires a zero-point. This leads to a hole in the atomic density, 
the width of which is only limited by the intensity ratios between the incident laser beams 
due to the geometric nature [25] of the process. The formation of the hole is accompanied 
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Figure 1: a) The level scheme for the three laser beams Oj (i = 1,2,3). b) The sequence 
of laser beams being swept through the BEC involves a preparation stage Q2 — » ^i and 
a final stage Qi — > Q 2 which engineers the phase and density of the BEC to produce a 
soliton. 

by a step-like (infinitely sharp) ir phase slip in the atomic wave-function when crossing 
the zero-point. The method is particularly useful for creating multicomponent (vector) 
solitons of the dark-bright form as well as the dark-dark combination. In addition it is 
possible to create in a controllable way two or more slowly moving dark solitons close to 
each other for studying their collisional properties. 



2 Formulation 

2.1 Outline of the proposed setup 

Consider a cigar-shape atomic BEC elongated in the z-direction. To create solitons in the 
BEC, we propose to sweep three incident laser beams across the condensate. The laser 
beams interact with the condensate atoms in a tripod configuration [25,26], i.e. the atoms 
are characterized by three ground states |1), |2), |3) and an excited state |0). The j-th 
laser drives resonantly the atomic transition between the ground state \j) and the excited 
state |0), see Fig. [Th. Initially the atoms forming the BEC are prepared in the hyperfine 
ground state |1). Subsequently the lasers are swept through the BEC in the x-direction, 
i.e. perpendicular to the longitudinal axis z of the condensate. 

The sweeping process is made of two stages depicted in Fig. [lb. In the first stage the 
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lasers 1 and 2 are applied in a counter-intuitive sequence to transfer adiabatically the atoms 
from the ground state |1) to another ground state |2). If an additional laser 3 is on during 
the first stage, a partial transfer of atoms from the ground states |1) to |2) is possible [26]. 
In that case a coherent superposition of states |1) and |2) is created after completing the 
first stage. In the second stage, the lasers 1, 2 and 3 are applied once again to transfer 
atoms from the state |2) back to the state |1) and from the state |1) to the state |2). If 
the amplitude of one of these lasers fii or fi 2 changes the sign at z — z , the BEC picks up 
a it phase shift at this point after the sweeping, and a soliton can be formed. This is the 
case e.g. if one of the beams is the first order Hermite- Gaussian beam centered at z — 0. 

It is important to realize that at least two laser fields are needed to complete the 
adiabatic transfer of population between the ground states. Therefore the adiabaticity can 
be violated in the vicinity of the point z = z where one of the Rabi frequencies Qi or Q 2 
goes to zero. Inclusion of the third (support) laser 3 helps to avoid such a violation of the 
adiabaticity. In fact the atoms would experience absorption in the vicinity of z = zq if the 
support laser 3 was missing during the second stage. 

It should be mentioned that there are similar previous proposals for creating vortices in 
a BEC via the two-laser Raman processes involving the transfer of an optical vortex to the 
atoms [27-29]. In these schemes the lasers are far detuned from the single-photon resonance 
to avoid the absorption at the vortex core. In our scheme the lasers are in an exact single- 
photon resonance, so the use of the third (support) laser is essential to avoid the losses. 
An advantage of the resonant scheme is that an efficient and complete population transfer 
is possible between the hyperfine ground states, whereas in the non-resonant case only a 
fraction of population can be transferred [29] . 

2.2 Hamiltonian for a tripod atom 

Let us now provide a quantitative description of our scheme. The j-th laser beam is 
characterised by the complex Rabi frequency Qj = Qj exp(?kj • r + iSj) , with j = 1, 2, 3, 
where flj is the real amplitude, the phase being comprised of the local phase kj r as well as 
the global (distance-independent) phase Sj. In what follows, the Rabi frequencies f2 2 and 
fi 3 are considered to be positive: f2 2 > 0, f2 3 > 0. Yet, the Rabi frequency f2 x is allowed 
to be negative. This makes it possible to include an additional tt phase shift in the spatial 
profile of the first beam when crossing the zero-point at z = zq. 

The electronic Hamiltonian of a tripod atom reads in the interaction representation: 



The tripod atoms have two degenerate dark states \D\) and \D 2 ) of zero eigen-energy 
(H e \D n ) = 0) containing no excited-state contribution [25,26], 



H e = -/i(fti|0)(l| + 2 |0><2j + Q 3 j0><3|) + H.c. 



(1) 




(2) 



(3) 
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where = \j) exp(i(k 3 — kj) ■ r + i(S 3 — Sj)) (with j = 1, 2, 3) are the modified atomic 
state- vectors accommodating the phases of the incident laser fields, ( = Q1/Q2 is the ratio 
between the Rabi frequencies of the first and second fields, and £j are the normalised Rabi 
frequencies (j = 1,2,3), 

is = % n = v^? + ^ + "I (4) 

with £3 > and —00 < C < +00. The atomic dark states |£)i) and \D<2) depend on the 
centre of mass coordinate r through the spatial dependence of the Rabi frequencies Qj and 
state-vectors 



2.3 General equations of motion 

The full atomic state- vector of a multicomponent BEC is |$(r, £)) = X^=i \j)^j{ r i t), where 
the constituent wave functions t) describe the translational motion of the BEC in the 
internal state \j) of the tripod scheme. The wave functions ^(r, t) obey a multicomponent 
Gross-Pitaevski equation of the form 



z^|$(r,t)> 



2M 



V 2 + H e + V 



(5) 



where H e from Eq. ([T]) describes the light-induced transitions between the different internal 
states of atoms. The diagonal operator 



v= £ (^ + ^l^r)b')0'l 

l>j=0 



(6) 



accommodates the trapping potential Vj(r) for the j-th internal state, as well as the non- 
linear interaction between the components j and I characterised by the strength gji = 
4irh 2 aji/m, with a,ji being the corresponding scattering length. 



3 Time-evolution of the atom-light system 
3.1 Adiabatic approximation for the dark states 

We shall apply the adiabatic approximation [25,30,31] under which atoms evolve within 
their dark-state manifold during the sweeping. This is legitimate if the total Rabi fre- 
quency fl is sufficiently large compared to the inverse sweeping duration r s ~g . The full 
atomic state-vector can then be expanded as: |$(r, t)) = J2n=i ( r > t)\D n (r, £)), where a 
composite wavefunction (r) describes the translational motion of an atom in the dark 
state \D n (r, t)). The atomic centre of mass motion is thus represented by a two-component 
wave-function 

* = (»?»)■ (7) 
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obeying the following equation of motion [25]: 



**** 



2M 



-ihV - A) 2 + V{r) + <f>-/3 



where the effective vector potential A and the matrix f3 are the 2x2 matrices appearing due 
to the spatial and temporal dependence of the dark states: A„ im = ih(D n (r, t)\ VD m (r, t)) 
and /3 nm = ih(D n (r,t)\d/dtD m (r,t)). The former A is known as the Mead-Berry con- 
nection [32,33], whereas the latter matrix j3 is responsible for the geometric phase [34]. 
The 2x2 matrix cf) is the effective trapping potential (explicitly presented in Ref. [25]) 
appearing due to the spatial dependence of the dark states. Assuming that all three beams 
co-propagate (ki ~ k 2 « k 3 ), the effective vector potential [25] reduces to 



(9) 




Lastly, the 2x2 matrix V originating from the operator V, Eq. ([6]), accommodates the 
trapping potential for the dark states [25] as well as the atom-atom coupling. 



3.2 Time-evolution during the sweeping 

Suppose the incident laser beams are swept through a trapped BEC along the x axis with 
a velocity v, as shown in Fig. [Tb. This can be done either by shifting in the transversal 
(x) direction the laser beams propagating along the y axis or by applying a set of laser 
pulses of the appropriate shape and sequence propagating in the x direction. In the latter 
case, the sweeping velocity v will coincide with the speed of light. In both cases the 
adiabatic dark states depend on time in the following way: \D n (r,t)) = \D n (r')), where 
r' = (x',y,z) = (x — vt,y,z) is the atomic coordinate in the frame of the moving laser 
fields. Let us assume that the time, r sweep = d/v, it takes to sweep the laser beams through 
a BEC of the width d, is small compared to the time associated with the BEC chemical 
potential t im = h/ji which is typically of the order of 10~ 5 s. In that case one can neglect 
the dynamics of the atomic centre of mass during the sweeping. Consequently the time 
evolution of the multicomponent wave-function during the sweeping is governed by the 
matrix-term (3 = —vA x featured in Eq. ([51) , giving 

ihdtV = vA x $> , (10) 

where the A x the effective vector potential along the sweeping direction. 

In passing we note that the subsequent time evolution of the BEC after the two-stage 
sweeping will be described by the general Gross-Pitaevski equation (jHJ) with the light fields 
off (H e = 0), as we shall do in Section HI 

Returning to Eq. ffTUl) . since vA x commutes with itself at different times, one can relate 
the wave-function *&(t) at a final time t — tf to the one at the initial time t — ti as 

*(r,t / )=exp(-ze)*(r,^), (11) 
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where the exponent G is a 2 x 2 Hermitian matrix 

9 = I A x {r - vt)vdt = I r A x (r')dx' . (12) 
a Jti ri Jxf 

and the integration is over the sweeping path r' = (x — vt,y,z) from Xf = x — vtf to 
Xi = x — vti. In most cases of interest the initial and final times can be considered to be 
sufficiently remote, so that the spatial integration can be from Xf — — oo to Xi = +oo. 

3.2.1 The first stage 

Let us now analyze the proposed two-stage setup in more details. In the first stage both 
Rabi frequencies fli and 2 are positive. The lasers 1 and 2 are applied in a counterintuitive 
order (see Fig. [Tb), where the ratio ( = f2i/f2 2 changes from = to C(^/) = +oo. 
On the other hand, the laser 3 is dominant for both the initial and final times where 
£3 = Q3/Q = 1. Initially the BEC has the wave-function \l/(r) and is in the internal 
atomic ground state |1) which coincides with the first dark state at the initial time t' iy i.e. 
\Di(r, t'j)) = |1). The full initial atomic state-vector is therefore |$(r,t-)) — ^( r )|-Di( r , t'i))- 
This provides the following initial condition for the multicomponent wave-function: 



= ( 7 ; j ■ (13) 

Equations ([9]) and ( fTTl) - (|T3l) yield the multicomponent wave-function after the first stage 
where 

/■+00 SarctanC , , 

is the mixing angle between the dark states acquired in the first stage. 

Suppose we have the following laser beams. The second beam Q2 is the Gaussian beam 
characterised by a waist a z in the z direction. The beam is centered at x' = x + A in the 
sweeping direction and at z = in the z direction, 

n 2 = Ae~ z2/a '- {x '-^ A)2/a '. (16) 

The first beam Q\ is characterised by the same amplitude A, the same waist a z and the 
same width a x . Yet it is centered at x' = x — A in the sweeping direction, where 2 A 
is the separation between the two beams. The beam waists should be of the order of the 
condensate length (or larger) in the ^-direction, so that the whole condensate is illuminated 
by the beams. 

The third beam is considered to change little along the sweeping direction x. Further- 
more it has the same width a z in the ^-direction as the first two beams 

fi 3 = AKe~ z2/a2 . (17) 
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Figure 2: Dependence of the mixing angle (3 on the relative amplitude of the third beam k. 
The spatial separation between the first and the second beams is taken to be 2 A = 1.2a x 
where a x is the width of the beams in the sweeping direction. 

The first stage is aimed at creating a superposition of states |1) and |2). Since we take 
all the beams to be the Gaussian beams characterized by the same widths a z , the Rabi 
frequency ratios f^/^i and ^3/^1 have no z-dependence. As a result the acquired mixing 
angle (3 has no ^-dependence, i.e. it is uniform along the BEC. The magnitude of j3 depends 
on the relative intensity of the third laser. If the third laser is weak (£ 3 = Q 3 /Q — > at the 
crossing point where ( = Q1/Q2 = 1), the mixing between the states |1) and |2) is small: 
/3< 1. On the other hand, if the Rabi frequency SI3 is comparable with Q\ and Q2 at the 
crossing point where ( = — 1 ; the mixing can be close to its maximum: (3 ~ 7r/4. 

In this way, one can control the mixing angle by changing the intensity of the third beam, 
as one can see from Fig. [2J 

3.2.2 The second stage 

In the second stage the Rabi frequency Qi can be both positive and negative depending on 
the transversal coordinate z. The laser 1 is now applied first, so that the ratio ( = Q1/Q2 
changes from ((ti) = ±00 to C(^/) = in the second stage. Again the third laser dominates 
for the initial and final times: Q3/Q = 1. The second stage takes place immediately after 
completing the first stage, so the multicomponent wave-function of the first stage 
serves as an initial condition for the second stage. 

Equations ffTT]) . ffT2"]) and (Till) together with and §5§ yield the total state vector after 
the second stage: 

|<3>(r,£/)) = |l)\I/(r) (sin7cos/? - e lI/12 cos7sin/?) 
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-|2)*(r) (cos 7 cos + sin 7 sin/?) . (18) 

where Vvi = Si — S2 + S' 2 — S[ is the phase mismatch between the Rabi frequencies Qi and 
Q2 in the first and second stages. The resulting mixing angle acquired in the second stage 
is 

d arctan £ 



7 = lz 



/" +0 °/, j. xC arctan C , . . . 

If the first and second lasers are weak — > 1 at the crossing point where C = ^1/^2 = 
1), the mixing angle is small 7 2 <C 1. On the other hand, if first and second lasers are 
strong at this point, we have 7 2 — > =F7r/2. The change in sign of 7 2 will introduce a phase 
shift which is needed to create solitons. 

In the second stage the first beam Qi is a first-order (in the z direction) Hermite- 
Gaussian beam centered at z = and x' = x + A 

Ql = At-J/oi-V-z-Wo*^ (20) 
B 

where z = ±5 represents a distance where Q\ = ±0,2 for x' = x. In most cases of interest 
the distance B is much smaller than the waist of the beams: B ^ a z . The second beam 
Q2 is the ordinary Gaussian beam centered at z = along the BEC and x' = x — A in the 
sweeping direction 

n 2 = Ae- z2/a *- ix '-* +K)2/a *, (21) 

where 2A is the separation between the two beams. The ratio between the Rabi frequencies 
reads then 

C = ^ = |^- fi)/oS , (22) 
Equation (122j) provides the following limiting cases: 



Finally let us determine the crossing point where Qi = ^2- Using Eq. (1221) . the condition 
|C| — 1 yields the crossing point x' = x' cr for a fixed z coordinate: 

a 2 z 

x' CI = x + ^\n — . (24) 
4A B v ' 

Specifically, if z = B, the crossing point is: x' CT = x. Since B <^ a, the Rabi frequencies at 
z = B and x' = x are: 

fix = n 2 w Ae~ A2/<T l (25) 

In the next subsection we shall analyse in more detail the multicomponent wave- function 
after completing the second stage. 
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Figure 3: Multicomponent wave-function after completing the second stage in the case 
where the second component is populated after the first stage {(3 = 0) and there is no 
phase mismatch between the lasers of the first and second stages (V12 = 0). The second 
and third laser beams are taken to be the Gaussian beams with equal widths o z . The first 
laser beam is the first order Hermite-Gaussian beam with same width a z . The parameters 
used are 2A/a x = 1.2, B/a z = 0.1 and k = 0.1. The wave function of the first (second) 
component is plotted in a solid (dashed) line. 



3.3 Multicomponent wave-function alter the sweeping 

Suppose that there is no phase mismatch between the lasers of the first and second stages: 
u 12 = 0. In that case Eq. (1T51) yields 

Mr,t f )) = *(r)[- sin( 72 - /3)|1) + cos( 7 , - /3)|2>] , (26) 

If (3 = 0, the second component is populated after the first stage. After the whole sweeping 
the state-vector then takes the form 

\${T,t f )) =*(r)[-sin 7 ,|l)+cos7,|2)]. (27) 

In this case the first component alters the sign at z = Zq where the Rabi frequency Qi 
or Q2 (and hence 7z ) crosses the zero-point. On the other hand, the second component 
is maximum at this point and symmetrically decays to zero away from this point. Such 
a multicomponent wave-function has a shape close to that of a soliton of the dark-bright 
form (see Fig. [3]) • This will indeed lead to the formation of such a soliton, as we shall from 
the analysis of the subsequent time-evolution presented in the next Section. 

On the other hand, (3 = 7r/4 corresponds to the case where both components are initially 
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Figure 4: Multicomponent wave-function after completing the second stage in the case 
where where both components are initially populated after the first stage {(3 = 7r/4) and 
there is no phase mismatch between the lasers of the first and second stages {1*12 = 0). The 
second and third laser beams are taken to be the Gaussian beams with equal widths o~ z . 
The first laser beam is the first order Hermite- Gaussian beam with same width a z . The 
parameters used are 2A/a x = 1.2, B/a z = 0.1 and k = 0.1. The wave function of the first 
(second) component is plotted in a solid (dashed) line. 

populated with equal probabilities. Thus we have after the sweeping: 

|$(r, t/)) = -*(r)[- sin( 72 - vr/4)|l) + sin( 7z + tt/4)|2)] . (28) 

In that case both components of the wave-function acquire a it phase shift in a vicinity of 
z = Zo where Qi = 0, as one can see clearly in the Fig. H] Note that the zero-points of 
each component are slightly shifted with respect to each other. This makes it possible to 
produce two component dark-dark solitons oscillating around each other, as we shall see 
in the following Section. 

If P = 7r/4 , yet there is a 7r/2 phase mismatch (vi 2 = 7r/2), Eq. ( TTBl reduces to 

|$(r 1 t / )> = -«r(r)e^i[|2>-i|l>] . (29) 

In that case both components are characterised by the same spatial modulation exp (ijz) 
and have a relative phase n/2 after the sweeping. Therefore both components initially 
have the same velocity distribution proportional to V^ z . Furthermore there is no hole in 
the atomic density of neither component after the sweeping, similar the case in the phase 
imprinting techniques. 
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In this way, the creation of solitons can be controlled by changing the mixing angle (3 
and the phase mismatch v 12 



4 Subsequent dynamics and soliton formation 

The optical preparation of the initial state of the two-component Bose-Einstein condensate 
described in the previous section, is fast compared to any characteristic dynamics in the 
Bose-Einstein condensate. This is the case if the time r sweep = d/v it takes to sweep the 
laser beams through a BEC of the width d, is small compared to the time associated with 
the BEC chemical potential r /t = h/fi which is typically of the order of 10~ 5 s. With the 
prepared initial state and for sufficiently low temperatures we can therefore describe the 
subsequent dynamics using a two-component Gross-Pitaevskii equation [12] 

ih-*! = [-—V 2 + V(z)+g n \y 1 \ 2 + g 12 \y 2 \ 2 ]y 1 (30) 

d h 2 

ih-q 2 = [-— V 2 + V(z)+g 22 \^ 2 \ 2 + g 12 \^ 1 \ 2 }^ 2 . (31) 

The external potential is here chosen to be quadratic in the ^-direction, 

V(z) = ^mu 2 z 2 , (32) 

where u is the trap frequency and m the atomic mass. The two-body interactions are 
described by 

9ii = i,j = {1,2} 33 

mi 

with the scattering lengths which represents the intra and inter collisional interactions 
between the atoms in the states 1 and 2. In Eq. f[3"3"j) we have introduced the effective 
cross-section S of the elongated cloud. Strictly speaking the elongated Bose-Einstein con- 
densate is three dimensional. If, however, the transversal trapping is sufficiently strong, 
the dynamics can be considered effectively one dimensional, as in Eqs. ( 1301) and ( I3T1) . This 
requires that the corresponding transversal ground state energy is much larger than the 
chemical potential of the condensate. We choose the normalisation as / cfe|\I/j(,2)| 2 = iVj, 
where iVj is the particle number in condensate % (i = 1,2). 

With the initial states from the previous section we can simulate the dynamics of the 
Bose-Einstein condensate. We consider a condensate with g\\ : g\ 2 : g 22 = 1.0 : 0.97 : 1.03 
where g\ 2 Ni = 286 and Ni = N 2 . The unit of length is J ^ and time is in units of u^ 1 . 
In figure we show the dark-bright soliton dynamics whose initial state is prepared by 
choosing /3 — and ui 2 = 0. The two-component system which has one dark soliton in 
component 1 and a bright soliton in component 2, is stable, i.e. the solitons are stationary. 
This shows that the initial state is indeed close to the exact soliton solution. If the initial 
state is prepared with f3 — 7r/4 and z/ 12 = 0, on the other hand, the dynamics is strikingly 
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Figure 5: The dark-bright soliton. The two figures show the one dimensional density as a 
function of time for component 1 and 2. The lighter (darker) colours correspond to higher 
(lower) atomic densities. 

different, see figureEJ In this case we create two dark solitons with opposite phase gradients, 
hence there is an oscillatory motion, sometimes referred to as a soliton molecule. Such a 
bound state is only stable if the soliton velocities are low [12] which is indeed the case here. 
Alternatively, with (3 = tt/A and z/ 12 = 7r/2, the solitons move in unison as shown in figure 
[7J The large oscillatory motion appearing in Fig. [7J stems from the fact that the condensate 
density is not homogeneous, hence the solitons experience an effective trap [11]. 

5 Conclusions 

In summary, we have proposed a new method of creating solitons in elongated Bose-Einstein 
Condensates (BECs) by sweeping three laser beams through the BEC. If one of the beams 
is the first order (TEM10) Hermite- Gaussian mode, its amplitude has a transversal it phase 
slip which will be transferred to the atoms thus creating a soliton. Using this method it is 
possible to circumvent the restriction set by the diffraction limit. The method allows one 
to create multicomponent (vector) solitons of the dark-bright form as well as the dark-dark 
combination. In addition it is possible to create in a controllable way two or more slowly 
moving dark solitons close to each other for studying the collisional properties. For this 
the first beam fli should represent a superposition of the zero and second order Hermite- 
Gaussian modes in the second stage. The soliton collisions will be considered in more 
details elsewhere. 
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Figure 6: The bound state dark-dark soliton. For sufficiently low initial soliton velocities 
the two dark solitons perform an oscillatory motion around each other. The figures show 
the one dimensional atomic density as a function of time for component 1 and 2. The 
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lighter (darker) colours correspond to higher (lower) atomic densities. 
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